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We study the effects of light scalar mesons on the skyrmion properties by constructing and ex¬ 
amining a mesonic model including pion, rho meson, and omega meson fields as well as two-quark 
and four-quark scalar meson fields. In our model, the physical scalar mesons are defined as mixing 
states of the two- and four-quark fields. We first omit the four-quark scalar meson field from the 
model and find that when there is no direct coupling between the two-quark scalar meson and the 
vector mesons, the soliton mass is smaller and the soliton size is larger for lighter scalar mesons; 
when direct coupling is switched on, as the coupling strength increases, the soliton becomes heavy, 
and the radius of the baryon number density becomes large, as the repulsive force arising from the uj 
meson becomes strong. We then include the four-quark scalar meson field in the model and find that 
mixing between the two-quark and four-quark components of the scalar meson fields also affects the 
properties of the soliton. When the two-quark component of the lighter scalar meson is increased, 
the soliton mass decreases and the soliton size increases. 

PACS numbers: 11.30.Rd, 12.39.Dc, 12.39.Fe, 14.40.Be 


I. INTRODUCTION 

More than 50 years ago, T. H. R. Skyrme proposed his 
pioneering idea that one may identify a certain nontrivial 
topological field configuration (soliton) in mesonic theory 
with the baryon [1]. Since then, this idea has been ex¬ 
tensively applied in particle physics, nuclear physics, and 
also condensed matter physics [2]. The originally pro¬ 
posed Skyrme model include only the degree of freedom 
of the pion. Motivated by the essential roles of meson res¬ 
onances such as the rho meson, omega meson, and light 
scalar mesons in nuclear physics, researchers extended 
the Skyrme model to Skyrme-type models to include such 
resonances [3, 4]. 

Recently, especially after the appearance of holo¬ 
graphic models of quantum chromodynamics (QCD), the 
effects of some resonances on the skyrmion properties 
were widely explored. Some key points related to the 
present work are summarized as follows: (i) by dimen¬ 
sionally deconstructing a Yang-Mills theory in a flat four¬ 
dimensional space in the large Nc limit, it was found that 
as more isovector vector resonances are included, the soli¬ 
ton becomes lighter [5]; (ii) by using a chiral effective 
theory of vector mesons based on the hidden local sym¬ 
metry (HLS) approach written up to the next-to-leading 
order, the omega meson that enters the chiral effective 
theory through the homogeneous Wess-Zumino term was 
found to provide a strong repulsive force [6, 7]; (iii) a 
scalar meson that is regarded as the Nambu-Goldstone 
boson associated with spontaneous breaking of the scale 
symmetry, i.e., the dilaton, was introduced to investigate 
its attractive effect [8-10]. It was found that the dilaton 
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can provide an attractive force of about 100 MeV in the 
model, which can reproduce chiral restoration in matter. 

In this paper, we explore the effect of scalar mesons 
made of two quarks and four quarks on the mass and 
size of the skyrmion by first constructing a chiral effective 
model including the pion, the rho meson, the omega me¬ 
son, and two-quark and four-quark scalar mesons. Then, 
determining the parameters of the model through meson 
dynamics, we study the effects of the masses and con¬ 
stituents of the scalar mesons, as well as the strength of 
an interaction among the vector and scalar mesons, on 
the mass and size of the skyrmion. 

Our findings can be summarized as follows: 

1. When we switch off mixing between the two-quark 
and four-quark states in the model, the four-quark 
scalar state decouples from the model. The light 
scalar meson is a pure two-quark state. In such a 
case, we find that: 

(a) The soliton mass and size decrease with in¬ 
creasing mass of the light scalar meson. 

(b) The soliton mass depends on the coupling 
strength between the light scalar meson and 
the vector mesons. The result shows that the 
soliton mass and the size of the baryon number 
density increase when the coupling strength 
increases. 

2. When we switch on mixing between the two-quark 
and four-quark states in the model, we find that, 
when the two-quark component of the lighter scalar 
meson is increased, the soliton mass becomes small 
and its size becomes large. 

This paper is organized as follows. The next section 
describes how we construct the chiral effective model that 
we use in this paper and explains its typical features in 
relation to this work. In sec. 3, we calculate the skyrmion 
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properties numerically and analyze the effects of scalar 
mesons on them. The last section discusses the results 
we obtained and presents our conclusions. In the ap¬ 
pendixes, we give some details of the calculation. 


II. THE MODEL 

We start by constructing a chiral effective model in¬ 
cluding two-quark and four-quark scalar states. Here, 
we limit our discussion to the mesons, including up and 
down quarks. In our model, we introduce a two-quark 
field M( 2 ) that is schematically written at the quark level 
as 


~ QRiqLj , ( 1 ) 

where i,j = 1,2 are the flavor indices. This, at the 
hadron level, is rewritten in terms of the isosinglet a field 
and isotriplet tt fields as 

M(^2) = ^{cr + in-T) , ( 2 ) 

where t are the Pauli matrices. Under chiral transfor¬ 
mation, the matrix M( 2 ) transforms as 

M(2) —t gLM(^2)9R , (3) 

where gL,R € SU(2)i^/{. Furthermore, we can impose 
Z 2 symmetry as a remnant of the U(l)^ transforma¬ 
tion, under which and transform as q^ —t —Ql 

I 

Uo = ATr (M( 2 )M 1 ' 2 )M( 2 )M 1 ' 2 )) 
-I- V 2 A ^det (Af( 2 )) -|- det 


We adopt the simplest form of the explicit chiral symme¬ 
try breaking potential Vsb: 

^SB = - Tr , (9) 

where x = and B and Ai are a constant with 

dimension one and the quark mass matrix, respectively. 
In this analysis, we consider an isospin-symmetric case, 
that is, Ai = diag(m, fh), where m = (m„ -I- md)l2,. We 
use the pion mass as an input to determine the relevant 
parameters, so Vsb is expressed as 

^SB = - luml Tr + M( 2 )) . (10) 

In this paper, we introduce the rho and omega mesons 
based on the HLS [12, 13]. For this purpose, we take the 
polar decomposition of M( 2 ) as M( 2 ) = where 

and in the unitary gauge of the HLS, are expressed 


and qR —> qR, respectively. Therefore, the hadron field 
M( 2 ) transforms as 


^(2) -t -M(2) . (4) 

We introduce a four-quark field that is schematically 
written at the quark level as 

(j} qLiqLjf"^ qRqR^ki ■ (5) 

In this work, we regard (j) as a real held by considering 
only the real part of Eq. (5). Under the Z 2 transforma¬ 
tion, (j) transforms as 


<j) ■ (6) 

On the basis of chiral symmetry and Z 2 symmetry, the 
Lagrangian is written in the form 

£ = Tr (a^M(2)5'^Mj2)) + 

- (Vb - Uq) - (UsB - Usb) , (7) 

where Vq represents the meson potential term, and Vsb 
denotes the explicit chiral symmetry breaking term, 
which has values of Vq and Usb, respectively, in vacuum. 
We note that, in this analysis, we include only the non¬ 
derivative interactions. We require that the number of 
helds in each vertex derived from the potential Vb is less 
than or equal to four, and the number of quarks included 
in each vertex is less than or equal to eight [II]. Then, 
the potential is expressed as 



by the pion fields as 

. ( 11 ) 

Then, we write the Lagrangian (7) as 
B -f cr^ Tr (a_L^a([) -h 

- (Vb - Uo) - (UsB - Usb) -f Cy , (12) 

where 

^ [d^^R • • 4) . (13) 

and the newly introduced term Cy represents the La¬ 
grangian for the vector meson component, which will be 
explicitly given later. The potentials Vb and Usb are 
rewritten as 

Uo = ^Acr^ - -h -h -^Aa‘^4>, 




3 


= - \mlUaTi- {U + U^) , (14) 

where U = 

The stationary conditions of the potential of the sys¬ 
tem are obtained by taking the first derivative of the 
potential = Vb -I- Vsb with respect to a and (j)\ 


dV 1,0 2 , 

= -Acr'^ — mia + V2Acr(f> 
oa 2 

9V 2^ 1.2 




(15) 

(16) 


After a suitable choice of the model parameters, cr ac¬ 
quires its expectation value cTvac, which is identified as the 
pion decay constant. Through the interaction between 
the two-quark and four-quark scalar states, a nonzero 
A leads to a nonzero expectation value of 4>, i.e., 4>va,c- 
The physical scalar meson fields are obtained from the 
fluctuations with respect to the relevant expectation val¬ 
ues through cr = Jtt + a and (f) = 4>vac + </>, where 
A. 

Taking the second-order derivative, one obtains the 
mass matrix for cr and (j) as 


2 777^ 

''V0 

rnl 


m- 

,2 


\f^ + m 

V2AU 


V2AU 


ml 


(17) 


- i[v^, K] , 


( 21 ) 


with 


V,= 


= — ( 
72 V 


Pt 




Here g is the gauge coupling constant of the HLS, Ohis 
and So are real dimensionless parameters, and F is a 
constant of dimension one. When cr has its vacuum ex¬ 
pectation value (VEV), a^ac = f■K^ the vector meson 
mass is expressed as my = ahis5^(so/^ + (1 ~ So)F‘^), 
where OhisS^so/^ represents the mass obtained by spon¬ 
taneous chiral symmetry breaking, and ahis5^(l ~ so)F'^ 
is the chiral invariant mass. In this sense, sq accounts 
for the magnitude of the vector meson mass coming from 
spontaneous chiral symmetry breaking. In this paper, we 
take F = so the vector meson mass is expressed as 
my = ahis5'^/7 which is the standard form given in the 
HLS [12, 13]. Note that for £vb, the term linear in the 
cr field is excluded by Z 2 symmetry, and sq can be both 
positive and negative. 

The intrinsic parity-odd terms of the Lagrangian Tanom 
in Eq. (19) read [13] 



(23) 


The physical states /500 and /1370 are the mixing states 
of cr and 0 through the rotation 


where M'^ represents the four-dimensional Minkowski 
space, and 


{ /500 \ { cos6» -Sind A f a\ 

\fi37oJ~\sme cosejy^j^ 

where 6 is the mixing angle. In Appendix A, we show 
some relations among the parameters and the physical 
quantities that are useful for this work. Erom Eq. (18) 
one can see that when cos d —>■ 0, the physical state /500 is 
dominantly a four-quark state and /1370 is dominantly a 
two-quark state, whereas when cosd —>■ 1 , /500 is almost 
a two-quark state, but /1370 is almost a four-quark state. 

Now, we are in a position to specify the vector meson 
component of the Lagrangian. Here we use the following 
form: 


Ay {aLfj.aLvaLaaRp — aRpaRuaR^aLp) , 
C-2 = Ay {aLpOtRvaLaOLRp) , 

A-s = Tr [Fy^i, {ollgOlrp — aR^aRp)] , (24) 

in terms of the I-form and 2 -form notation with ur = 
(A|| ~ Q;_l, aR = d|| -|- a_L, and Fy = dV — iV^. The 
low-energy constants ci,C 2 , and C 3 , which are difficult 
to fix from fundamental QCD, are usually estimated 
phenomenologically with respect to the experimental 
data [13]. In this work, we focus on the effect of scalar 
mesons on the skyrmion properties; therefore, we choose 
the parameters ci = — C 2 = — 2/3 and C 3 = 0, which 
provides [14]. 


Cv = Avo + -Canom , (19) 

where Avq is the Lagrangian with only intrinsic parity- 
even terms, whereas £anom contains the intrinsic parity- 
odd terms relating to the V{2)r x U(2)7{ chiral anomaly. 

Explicitly, the intrinsic parity-even Lagrangian Avo is 
expressed as 

Avo = ahis(soCT^ + (1 - so)F'^) Tr(d||77d{|‘) 

-^AY{Vp,Vn, ( 20 ) 

where 0 : 11 “ and Vp^, are defined as 

^ ■ Cl) ~ j 


III. NUMERICAL RESULTS FOR THE 
SKYRMION 

In this section, we make a numerical study of the ef¬ 
fects of scalar mesons on the soliton mass and also the 
root-mean-square (RMS) radii of the baryon number den¬ 
sity and energy density. In Table. I, we summarize the 
parameters of the model used in this paper. 

A. The ansatz 

To study the properties of the soliton obtained from 
the Lagrangian (12), we take the standard parametriza- 
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u 

92.4 MeV 

ruTT 

139.57 MeV 

N, 

3 

Cl -1- C2 

0 

Cl - C2 

-4/3 

C3 

0 

9 

5.80 ±0.91 

Q-hls 

2.07 ±0.33 


TABLE I: Values of the model parameters. 


tions for the soliton configurations. Following Refs. [1, 
14], we take the ansatz for the tt, p and w helds as 

U = , (25a) 

_ ^ ^ (25b) 

gr 

uj^ = tV(r)5^o • (25c) 

We parametrize the scalar meson fields a and ip as 

cr =/tt (1 + d(r)) , (26a) 

(p = (pvac (l + (p{r)) , (26b) 

where (t(r’) and 0(r) are dimensionless functions. Sub¬ 
stituting Eqs. (25) and (26) into the Lagrangian (12), one 
obtains the equations of motion for the profile functions 
F{r), G{r), W{r), a{r), and ^(r). The detailed expres¬ 
sions are given in Appendix B. For the solutions with the 


baryon number B = 1, the wave functions E(r), G{r), 
and W{r) satisfy the following boundary conditions: 

E(0) = TT, F(oo) = 0, 

G(0) = - 2, G(oo) = 0, 

W'{Q) = 0, IV(oo) = 0. (27) 

The boundary conditions for the profile functions of the 
scalar mesons are determined as follows: (I) as r —oo, 
the VEVs of the scalar meson fields should be repro¬ 
duced, so both a{r) and (p{r) must vanish, i.e., 

(t(oo) = 0 , 0 (oo) = 0 ; (28) 

(II) the second-order derivatives a” and (p” should be 
nonsingular as r —>■ 0, in view of the last term in Eqs. (B 6 ) 
and (B7); equivalently, 

ct'(0)=0, (^'(0) = 0. (29) 

To study the RMS radius of the baryon number den¬ 
sity of the system, we specify the baryon number current 
by introducing the external gauge field Vfi of the 1 /( 1 ) 
baryon number. Then, the baryon number current is 
obtained by taking a functional derivative of the total 
Lagrangian (12) with respect to V^. After an explicit 
calculation, the baryon number density Bq, which is the 
zeroth component of the baryon number current, is ex¬ 
pressed as 


Bo = - {f^g'^r'^ahisW [socr^ -|- 2soa’ + l] 

+ F' \a 2 — 2G (—a 2 -I- 0:3 -I- 02 cos F) + 02 cos^ F — 2 q !2 cos F + (02 — 03 ) G^] 
-2a3sinFG'-haisin2FF'}--^F', (30) 

where ai, 0 : 2 , and 03 relate to ci, C 2 , and C 3 , respectively, through Eq. (B2). By using the equation of motion for W 
in Eq. (B5), the baryon number density is simplified as 

Id / 4 a 3 sinF(cosF — G — 1) 2r'^W' sin(2F) F \ 

° dr \ ig 3g Stt^ 47 r^y ’ 


which agrees with that obtained in Ref. [15]. Note that, 
by combining Eqs. (31) and (27), one can show that the 
baryon charge is correctly normalized as dr^Bo(r) = 

1 . 

In this analyse, we take ci -I- C 2 = 0, ci — C 2 = 3/4, and 
C 3 = 0; therefore, Eq. (30) is reduced to 

Bo= - [sqct^ -f 2 soct -t l] . (32) 

In our numerical calculation, in addition to the soli¬ 
ton mass, we also consider the RMS radius of the baryon 


number density, , and the RMS radius of the en¬ 

ergy density, {r'^)^^, which are defined as follows: 

^ j d^rr'^Boir) , 

{r^)E^ = d3rr2Msoi(r), (33) 

where Mgoi is the soliton mass, and Msoi{r) is the corre¬ 
sponding energy density given in Eq. (Bl). 
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B. Effects of scalar mesons 

Now we are in a position to show how the scalar mesons 
affect the skyrmion properties. 

First, we consider A = 0. In this case, the four- 
quark field (j) decouples from the other mesons. The La- 
grangian (12) is reduced to 

= ^dfj,ad^a + Tr(a_L^a([) 

- (14 - K) - (4sb - Fsb) + Cv , (34) 

where 14 = The parameters A and 

are related to the masses of cr and tt through 

2 2 I 3 , .2 

-m2 + -XU , 

ml= -ml + ■ (35) 

Here we take the empirical value of m,r as an input and 
regard TOo- as a parameter of the model. In addition, 
So in is an unfixed parameter. In the following, we 
study the relation between these two parameters and the 
skyrmion properties. 

a. The effects of rria on the skyrmion properties. 
We start our discussion from the standard HLS La- 
grangian [13], which is obtained by taking sq = 0 in this 
analysis. 

We show the d epende nce of the soliton mass msoi 
and the radii U (r ^)b and \/(U)e in Fig. 1 and the pro¬ 
file functions in Fig. 2. From Fig. 1 we find that, in¬ 
creasing ttIct, the soliton mass ttIsoI increases. This ten¬ 
dency can be understood from the profile function of a 
in Fig. 2; when the a is heavy, the magnitude of a is 
small. Because ct < 0, as mo- increases, the magnitude of 
a = + O') becomes large; consequently, the contri¬ 

bution of cr^ Tr(a_L^a([() to the skyrmion mass becomes 
large. Physically, for a heavy cr, the attractive force pro¬ 
vided by the cr is small; as a result, the soliton mass is 
large. 

In a ddition . Fig. 1 tells us that when Wcr is large, the 
radii ■>/(r ^)e and a/ (r ^)b are small. This result can be 
understood from the profile functions plotted in Fig. 2; 
when mcr is large, the profile for a becomes narrow. Gen¬ 
erally, the cr supplies the attractive force, whereas the u! 
meson supplies the repulsive force. A narrow a needs a 
narrow w to balance the attractive and repulsive forces. 
The shapes of a and ui dominate the shape of the energy 
distribution. As a result, the narrow cr and uj make the 
radius a/( r ^)^ small. As for the radius of the baryon 
number density, from expression Eq. (32) with sq = 0, 
one can easily see that a narrow W gives a small radius 
\/ B- To summarize, when mo- is large, the effective 
range of the attractive and repulsive forces supplied by 
the a and uj, respectively, are small; i.e., the soliton size 
is small. 

b. The effects of sq on the skyrmion properties. 
The parameter sq controls the chiral invariant vector me¬ 
son mass. Here, we show how sq influences the properties 



FIG. 1: mo dependence of the soliton mass and radii 
with A = Sq = 0. 



r (fm) 


FIG. 2: mo dependence of the profile functions with 
A = So = 0 for mo = 1 GeV (green solid line) and 
m„ = 1.4GeV (red dash-dotted line). 


of the skyrmion by taking a typical value of the scalar 
meson mass, mo = 1.37 GeV. Our numerical results are 
plotted in Fig. 3. 

From Fig. 3, we see that, when the value of sq is in¬ 
creased, mgoi and U (r ^)b become large. According to 
the equations of motion for the vector mesons in Eq. (B5), 
OhisnS (/^(l — So) -I- socr^(r)) acts as the effective mass of 
the vector mesons inside the soliton. Because a{r) < f^^, 
as we stated above, the effective mass is smaller than the 
meson mass m^^ = ffaffffgujf-n in vacuum. This has two 
effects: on the one hand, the effective strength of the re¬ 
pulsive force mediated by the uj meson is strong for large 
So, which produces a large soliton mass, as the vector 
meson dominates the energy; on the other hand, the ef¬ 
fective range of the repulsive force supplied by uj is long 
for large so, which makes the RMS radius of the baryon 
number density U (c^)_b large, as the uj meson acts as 
the gauge boson of the 17(1)v baryon number symme- 
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try. The relation between sq and (r ^)b needs a more 
precise explanation; because the baryon number density 
is proportional to the effect of the uj meson, as can be 
seen in Eq. (32), the RMS radius of the baryon number 
density y/(r ^)b is large for large sq. 

The situation for the RMS radius of the energy density, 
E, is quite complicated. In this analysis, the nu¬ 
merical error is abou t 1-3%, so the third graph in Fig. 4 
implies that y/ (r ^)e is rather stable against changes in 
So- This implies that the contribution of w has the op¬ 
posite direction to that of tr, so they cancel each other. 


1.4 


^ '-35 

a 1.3 

. 


^* 


-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 


So 
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-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 
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FIG. 3: So dependence of the skyrmion mass and radii 
for nicr = 1.37 GeV. 



FIG. 4: Profile functions with So = — 0.5 (green line) 
and So = 0 (red dotted line). 


Next, we switch on the mixing between the two scalar 
mesons by setting A ^ 0 and then study how the mixing 
angle 9 affects the skyrmion properties. In our calcu¬ 
lation, we take two fixed values of so, i.e., So = 0 and 
So = — 0.5, to see the dependence of the mass and 
radii of the skyrmion on the mixing angle 9 defined in 


Eq. (18). Our results are shown in Figs. 5 and 6 . Note 



FIG. 5: Dependence of skyrmion properties on the 
mixing angle between two scalar mesons for so = 0 . 



FIG. 6 : Dependence of skyrmion properties on the 
mixing angle between two scalar mesons for sq = — 0.5. 


that we cannot find any stable solution in the regions of 
cos 9 > 0.6 for Sq = 0 and cos 9 > 0.25 for sq = — 0.5. In 
other words, the maximal value of cos 9 is small for small 
Sq. This could be understood as follows. The definition 
of the mixing angle in Eq. (18) implies that, for a large 
value of cos 9, /500 includes a large amount of the two- 
quark component a, which generates a strong attractive 
force. Because the total repulsive force does not change, 
the attractive force can exceed the total repulsive force, 
and the soliton collapses. A large value of cos0 corre¬ 
sponds to a strong attractive force; therefore, the soliton 
collapses at a certain value of cos 9. On the other hand, a 
decrease in sq means that the repulsive force mediated by 
the oj meson becomes weak, as we stated above. Because 
the soliton is stabilized by the balance between the at- 
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tractive and repulsive forces, the maximal value of cos 6 
is small for small sq. 

Note that the decreasing and increasing tendencies in 
Figs. 5 and 6 are the same, and they are actually consis¬ 
tent with those in Fig. 1. As we just stated, increasing 
cos 6 increases the magnitude of the two-quark compo¬ 
nent included in / 500 . Thus, the large cos 6 in Figs. 5 
and 6 corresponds to the small mass in Fig. 1. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this work, we first investigated the effect of a two- 
quark scalar meson on the skyrmion properties when mix¬ 
ing between the two-quark and four-quark states of scalar 
mesons is switched off. We found that as the mass of the 
two-quark scalar meson decreases, the soliton mass de¬ 
creases and the RMS radii of the baryon number density, 
5 , and the energy density, \/(r^)increase. We 
next studied the effect on the skyrmion properties of cou¬ 
pling between the scalar meson and vector mesons. It is 
found that as coupling becomes larger, the soliton mass 
and radii increase. We finally explored the skyrmion 
properties by considering the mixing structure of the two 
scalar mesons. We showed that, when the percentage of 
the two-quark component in the lighter scalar meson is 
increase d, the soliton mass becomes small, and the RMS 
radii (\/(r^)B and \/(r^)_E) become large. 

Our numerical results show that, even though the two- 
quark and four-quark scalar mesons are both included in 
the model, the soliton mass obtained here is still about 


300 MeV larger than that needed to reproduce the em¬ 
pirical values of the nucleon mass and A-V mass split¬ 
ting. Thus, to reproduce the baryon masses in the real 
world, the present model should be extended. There are 
several possible extensions, for example, inclusion of the 
effect of the dilaton as the Nambu-Goldstone boson as¬ 
sociated with the scale symmetry breaking of QCD; the 
effect of the terms of the HLS, which account for 

the TT-p interaction; and/or the effect of all the homo¬ 
geneous Wess-Zumino terms of the HLS, which increase 
the number of ir-p and/or TT-p-uj terms. We leave these 
possibilities as future projects. In addition, even though 
the present model cannot reproduce reality, it still merits 
application to dense matter systems to study the effect 
of the scalar mesons on the phase structure of nuclear 
matter. We will report this progress elsewhere. 
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Appendix A: The mixing matrix 

The masses of /500 and /1370 are 


^hoo = - I - 


Ul3 


\/ 2/^ (4A2 - Xml + + (”^4 - 1 ^ 1 ) + f^X + ml + , 

i Y^2/2 (4A2 - Xml + Xml) + + {ml - ml) ^ + f^X + ml + ml^ . 


(Al) 

(A2) 


The mixing angle 0 can be obtained as 


0 = arctan < 




/500 


m2 — m 


= arccos ^ 


/500 


2 2 

mi — mi 

/1370 /500 


From Eqs. (15), (16), (Al) and (A2), one gets 

Ki 370 -^4)("i4-™Lo) 


= 




x = 


m 


/soo 


2 2 2 

mi — mi — ml 

/l370 4 TT 


mo = 


(A3) with < na\ < rrij^ 



1 / m\_rin?i 


- 


(A4) 



















Appendix B: Soliton mass and the equations of motion for the profile functions F{r),G{r),W{r),a{r) and (j>{r) 


Substituting the ansatz defined in Eqs. (25) and (26) to the effective Lagrangian (12), we obtain the following 
expression for the soliton mass: 


pOO 

Msoi = Itt / drr^Msoi(r) 

Jo 

= -inj dr|ir^|4/^g^soahisVb^d(d-+ 2) - 4/^ (d-+ 1)^F'- 


+/^ 


,2 (ml - - ml) (d + 1)" + l) 


mt 


+ - 


(r 


2 (ml ml 

‘ /500 /1370 


- (™/500 + "^/1370 - - ml) (a + 1) 

3mlml) (d + 1)^ (ml - (ml - + l)^ 


+ 


(ml - ml^„„) (ml - ml^3^„) 


2-7 


4 

ml 


ml 

fl ( 


+ 8ml (tr + 1) cos{F) 


^%oo^h37o + 3"4lml 

ml 


■4W‘ 


/2 


} 


^|8/^5^ahisGsin^ (^y^ (sod^ + 2sod + l) - Aflg^a"^ sin^ (^^^ ((soUhis - 1) cos(F) - so-a-his - 1) 

-8/lg^dsin^ (^^^ ((soOhis - 1) cos(F) - soUhis - 1) + 2/l5^ahisG^ (sod^ + 2sod- + l) 

+8/lg^ahis sin^ (f) “ cos(2E) + flg^ + 2G'^| 

-03 ( 2 G {WF' - sin(E)lE') + G^IEF' + 2 sin(E) ((cos(F) - l)W' + WG')) 
+a2WF'{-cos{F)+G+lf FaiWF'sm^{F)- ^ 


2g^r 




with 


35^c , ^ 5^c , , ^ 

Q-l = , ^ T (<^l — -32) , 0-2 = 77~5'(-^i + <^ 2 ) , 0-3 = TTTJS^'i ■ 


16^2 V- 10^2 

The equations of motion for F(r), G(7’), IT(r-), d(r) and (j){r) are expressed as 

F” = ^ 


(Bl) 

(B2) 


/2r-2(l + d)2 

X | 2 G (/Isoahisd^ sinF + 2 /^soahisd sin F + /^Ohis sinF 

- 02 cos FIT' - 03 cos FIT' + a2VTG' - oglTG' + 02IT' - oaVT') 

- fla sin F (dsoUhis cos F — dsoOhis — 4 cos F — r^m^) — 2 /ld^ sin F (soOhis cos F — So-^his — cos F) 

- 2 flr'^F'a' - 2 flra‘^F' - 2 flraF' {ra' + 2 ) + 2 /lahis sinF - /^ahis sin( 2 F) 

- 2 /^rF' + flmlr^ sinF + /I sin( 2 F) — 202ITcosFG' + 203ITcosFG' 

- ioi cos( 2 F)lT' - 2o2 cos FIT' + ^02 cos( 2 F)VT' - 203 cos FIT' + 203 cos( 2 F)lT' + 202ITG' 

- 203ITG' + (02 - 03) G^IT' + ioiVT' + ^02lT'|, 

G" = - flg^ams (sod^ + 2 sod + l) (cosF - G - 1 ) - a^g^ (lTF'(cosF - G - 1 ) + 2 sin FIT') 

G (G 2 + 3 G + 2 ) 


+ a2g^WF\c,os F — G — 1) 

IT" = flg'^amsW (sod^ + 2 sod + l) - 


(B3) 

(B4) 


21T' 


+ ^{ - 03 (F' (2G(cosF + 1) + 2(cosF - cos(2F)) + G^) + 4sinFG') 
+ 02 F'(— cosF + G + 1)^ + 2oiF' sin^ (y) + 4)| > 


(B5) 
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^ (-S^soahisH^^ + + F'^ - ml) + {m%^^ + -mj- ml) 

+ 3 (™/ 5 oo + ™/i 37 o - ™4 - ^l) + 2 [ - g'^soausW^ + F'^ - m^(cos F - 1)] 

- {ml - m%J {ml^^^ - ml) {a + 1) 

(ct + 1) (-SsoOhisGsin^ {^) - SsoOhis sin'^ (-|) - 2soahisG^ + cos(2F) - l) 2a' 

^^2 j, ’ 

= ml {<f> - a {a + 2 ))- — . 


(B6) 

(B7) 
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